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ABSTRACT 

^ , We generalize the superposition principle for time-symmetric initial data of 

^ I black hole spacetimes to (anti-)de Sitter cosmologies in terms of an eigenvalue 



problem Agcf) = ^{Rg — 2A)0 for a conformal scale applied to a metric gij with 
constant three-curvature Rg. Here, Rg = 0,2 in the Brill-Lindquist and, respec- 
tively, Misner construction of multihole solutions for A = 0. For de Sitter and 
anti-de Sitter cosmologies, we express the result for i?^ = in incomplete ellip- 
o'l tic functions. The topology of a black hole in de Sitter space can be extended 

into an infinite tower of universes, across the turning points at the black hole 
{^jqI and cosmological event horizons. Superposition introduces binary black holes 

for small separations and binary universes for separations large relative to the 
^ _ cosmological event horizon. The evolution of the metric can be described by a 

Tij- I hyperbolic system of equations with curvature-driven lapse function, of alternat- 



(N 



o 



X 



ing sign at successive cosmologies. The computational problem of interacting 



O I black hole-universes is conceivably of interest to early cosmology when A was 

I large and black holes were of mass < |A~^/^, here facilitated by a metric which 

Q I is singularity-free and smooth everywhere on real coordinate space. 



Introduction 



The naultiho le solutions of Brill-Lindquist and Misner ( iMisnerl Il963t iBrill fc Lindquist 



1964J : ICookll200l[ ) arise out of a superposition principle in the Hamiltonian constraint of the 
Einstein equations for time-symmetric initial data. Here, the Schwarzschild line-element in 
coordinates (r, 6, (p) is transformed into an isotropic line-element described by a conformal 
factor 0(r), giving rise to a conformally flat spacetime in vacuum. This construction is 
remarkable, in allowing for multihole solutions with different extended topologies. To a 
black hole binary, the Brill-Lindquist solution attributes a three-sheet topology, whereas the 
Misner solution attributes a two-sheet topology to the same. 

The embedding of black holes in (anti-)de Sitter space is well-known for a single black 
hole with one-sheet topology. De Sitter space is of interest in view of considerable observa- 
tional evidence for a small but distinctly positive cosmological constant in the ACDM model 
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( Ide Bernardis et al.ll2000l : iHanany et al.ll2000l) with an expected improved uncertainty by the 
recently launched Planck satellite JeSA-SCI 2005 ). The stability and thermodynamics of 
the cosmological event horizon in the black hole-de S itter space is relevant to its potential 
role to early cosmology (jPaviesI 119871 : 1 Chambers! 119971 ) . 



Here, we consider the problem of extending space beyond the cosmological event horizon. 
We approach this problem in terms of the turning points at extrema of the circumference 
in the one black hole-de Sitter space. The result opens the possibility for novel topologies 
on scales that reach beyond the visible universe. Space of finite volume further facilites 
representation in a metric that is smooth everywhere. This is of direct interest to introducing 
spectral methods for calculating wave-templates for gravitational-wave observatories LIGO 
and Virgo as they are gradually improving their sensitivity. 

To begin, we recal the Brill-Lin dquist line-element of a single Schwarzschild black hole 
(IBrill fc Lindquistlll964j : ICookll200ll ). The event horizon of a black hole of mass E represents 
a turning point in the embedding in two asymptotically fiat sheets. 



over the donut ds\) 



- tanh2(A/2)rft2 + 4^2 cosh^(A/2)rf4, (1) 
l — x'^)d(p'^ with —00 < A < 00. Here, we transformed the 



dx 



+ 



Brill-Lindquist line-element with conformal factor 1 + ^ in spherical coordinates {p,9,ip), 



X = cos ^ by p = ^Ee^. The Mobius invariance A ^ —A comes with opposite signs in the 
lapse fuction = tanh(A/2) on either sheet, wherein the horizon surface corresponds to the 
simple zero = 0. 

By Liouville's theorem, spacetime singularities are inherent to any black hole spacetime 
which is asymptotically regular. They can be moved away into the complex plane as ([I]) 
illustrates by mapping p = —\E to A = 7ri (mod 27r), whereby the metric is analytic 
everywhere for all real and finite A, wherein A = ±00 represent coordinate singularities 
associated with asymptotic infinity on each sheet. 

The black hole singularity is not directly accessible by observation, whether in real or 
complex coordinate space. Even an observer in free fall onto a black hole never reaches the 
central singularity when considered in real coordinates, as in the Schwarzschild line element. 
Upon approaching the event horizon, time-at-infinity, t, becomes arbitrarily large, which 
signals evaporation of the black hole by Hawking radiation. The observer hereby traces a 
shrinking event horizon, never to penetrate it, during which time the black hole singularity 
diminishes in strength. Note that this result only uses evaporation of a black hole in a 
finite time-at-infinity, wherein the velocity of the black hole horizon is representative for 
the luminosity regardless of the details of Hawking radiation. Hawking radiation hereby 
introduces invariance to the cosmic censorship conjecture, by treating singularities in real 
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and complex coordinate space on equal footing. The same arguments shows that ([T]) defines 
a non-traversable wormhole. 

In de Sitter space, coordinate singularities associated with the infinite extend of space 
are avoided. The result for multihole solutions in de Sitter space, then, is a metric which 
is everywhere nonsingular and finite on the coordinate cover of the visible universe with no 
need for compactification. The cosmological event horizon represents an additional length 
scale, which can be of interest in considering black hole binaries with large separations. It 
introduces an additional turning point which opens a window for novel large scale extensions 
beyond. 

In §2, we derive the formalism for constructing multihole initial data in (anti-)de Sitter 
space. In §3 we give some illustrative examples for binary black holes and binary cosmolo- 
gies. In §4, we propose a hyperbolic s ystem of equations for their evo lution based on the 
3+1 Hamiltonian equations of motion (lArnowitt. Deser fc Misnerlll962l), where hyp e rbolic - 



ity generally facilitate s stabl e numerical implementation 



Nagy. Ortiz, fc Reulal (l2004j ): ICalabrese. Hinder fc Husal ( l2006l ) and references therein) by 



eg^ 



van Putten fc EardlevI ( 119961 ): 



ensuring a real dispersion relation and hence s tability whenever the Courant-Friedrichs-Lewy 
condition ( ICourant. Friedrichs fc Lewylll967| ) is satisfied. An outlook is included in §5. 



2. A superposition principle in (anti-)de Sitter space 

The line-element of a Schwarzschild black hole of mass m in a de Sitter space with 
cosmological constant A is 

ds^ = -(1 - 2m/r - -Ar^)dt^ + ^ . + r^dO^ + sin^ ed<f^ (2) 

3 1 — 2m/r — ^Ar^ 

When Am^ < | , a black hole of mass m can exist within the larger cosmological event horizon 
defined by A with a black hole horizon at coordinate radius r = ri and an event horizon at 
r = r2 > ri corresponding to the roots of the redshift factor = 0, A^^ = 1 — 2m /r — |Ar^, 
in ([2]). The roots are shown in Fig. ([1]) with the bifurfaction point at Am^ = |, including 
the non-physical real root < 0. The line-element ([2]) can be transformed to an isotropic 
line-element 

ds^ = -N^de + (P\dp' + p'de^ + p2 sin2 ed^'') (3) 
according to the conditions 

dv 

± I = = <l>'dp, r = ct?p. (4) 

1 — 2m/ r — ^Ar^ 



-4- 




Fig. 1. — Shown are the three real roots of the redshift factor in the Schwarzschild-de Sitter space 
as a function of Km? > 0. They define the coordinate location of the horizon of the black hole 
(ri) and of the cosmological event horizon (^2). The roots bifurcate into a pair of real roots at 
iW = 1/9. The maximum mass of a black hole in a de Sitter space is hereby l/SA-^/^ (jPodolskvl 



I999I ). Shown is further the numerical root rs < 0. 
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For A > the three real roots < < ri < r2 of N satisfy ri = 2m(l+e), r2,3 = ± y |-(1=f5) 
where e = |Am^ + 0(e^) denotes the horizon surface of the black hole and 6 = mJ ^ + 0{5'^) 



The solution to (jl]) can be expressed following a logarithmic transformation d\ = dp/ p. 
For A > 0, the Mobius transformation t = 2m /r (with no reference to time intended) gives 

\=( r - t\ ^ ^ F(7r/2,mi)-F(0,mi) 



in terms of the three roots = 2m/rj (^3 < < t2 < ^i), ""^i = f^z^; sin0 = y , 
wher e mi) denotes the incomplete elliptic function of the first kind (lAbramowitz fc Stegun 
19681 ). For A < 0, the only real root is r2- We proceed with r — ^"^ 



1-12?;' 



A = / (6) 
' v^4P3 - P/12 + 1/216 + AmVl2 



giving rise to the Weierstrass elliptic function v{\) = P{\, 1/12, —1/216— Am^/12) (lAbramowitz fc Stegun 



19681). 



In the resulting conformally flat approach, the resulting extension of de Sitter black hole 
spacetime is shown in Fi gs. (12HS1)- We note that the singular limit A = reduces to the 
familiar expressions (e.g. [Abrahams &: Pried (Il996[ )) 

A = In ^ — [r — m + Vr"^ — 2mr]^ , r = p ^1 + — ^ , p = - (i/r + Vr — 2m) ^ . (7) 
We interpret the result as follows. 

Theorem. The superposition principle for time- symmetric data in the Hamiltonian en- 
ergy constraint generalizes to (anti-)de Sitter space for eigenfunctions of the Laplace operator 
of a scaled metric with constant 3-curvature. 

Proof. We recall the conformal decomposition of the Ricci tensor for hij = (p'^Qij, 
whereby ^^^Rijih) = ^^'^Rij{g) + with dj = -20-i[A^i0 + gij^<P] + 2(p~^[3Di(j)Dj^ - 
gijD^(j)Dp(f)]. Following a reduction of the Ricci tensor from four to three dimensions, the 
Einstein equations for time-symmetric data give 

(^)R{h) = (3)i?(^) - 8(p-'Ag(P = 2A. (8) 

If the scaled metric gij produces a constant curvature Rg, the are eigenfunctions of the 
associated Laplace operator, 

Ag<p=l{Rg-2A)<f). (9) 
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Fig. 2. — For a single black hole in de Sitter space, the scale (p is periodic in the radial coordinate 
— oo < A < oo, representing a "tower" of cosmologies by successive continuations at the turning 
points defined by the black hole and cosmological event horizons. The results shown are computed 

1 /2 

for m = 1 and A = 0.001, and may be contrasted with </>=(! — e~^) in singular limit A = 0. 
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Fig. 3. — The extension of black hole-de Sitter cosmology of a single black hole-de Sitter spacetime 
can be realized by joining the eigenfunction solutions at successive turning points (solid-to-dashed 
transitions), representing black hole and cosmological event horizons (Hi,Ei). The results shown 
are computed for m = 1 and A = 0.001. 
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We can construct solutions by superposition 

= (10) 

/ij > 0, SyUj = 1, of different eigenfunctions 0j to the same eigenvalue, e.g., those translated 
in any one of the homogeneous directions. The Brill-Lindquist case corresponds to the fiat 
metric Qij = 6ij with Rg = and three homogeneous directions, whereas the Misner case 
corresponds to the donut metric of ([1]) with Rg = 2 and one homogeneous direction, both 
with A = 0. In case of the former, we have the coordinate density 



0j = ^, P= \J - V\? + {^2 - P2f + {^■i - P3)^ (11) 

associated with black holes at p (and r = in the corresponding Schwarzschild coordinates), 
where the singular limit A — )■ recovers 0=1 + ^. □. 

The eigenvalue problem (ITO!) reduces to a Helmholtz equation when is flat. It then 
suggests an association to a boson field with frequency uj = when A > 0. However, 

this association does not reflect scale invariance in the Riemann tensor, in that it depends 
only on the log of (p. 

A cosmological event horizon at finite distance alters the spectrum of eigenfunctions. 
Whereas the eigenfunctions r"^''^^^ Pi{x) to A0 = give convergence and divergence at 
asymptotic infinity on the two sheet embedding of the Schwarzschild black hole in the Brill- 
Lindquist case A = 0, the eigenfunctions of A0 = —\A(j) for A > are periodic across the 
black hole and cosmological event horizons in Figs. [2H3] and are hereby necessarily bounded 
and regular everywhere. 

A cosmological event horizon also alters the Hilbert space of radi ation states at large dis 



tances from a black ho le, and hence the details of Hawking radiation (IKanti. Grain fc Barrau 



20051 : IZeng et al.ll2008l ). This is particularly pertinent during inflation, when A is large. The 
evaporation time of black holes should remain finite, however, whereby the strong form of the 
cosmic censorship conjecture described in the introduction should continue to hold whenever 

A construction similar to (IMS]) may be pursued for the Misner two-sheet embedding 
of two black holes, based on eigenfunctions Ajjcp = 1(1 — a^A)0, where D refers to the 
donut line-element in ([1]) and a refers to the Misner length scale. A detailed consideration, 
however, falls outside of the present discussion. It ma y also be generized to extra d imensions 



on the basis of the Schwarzschild-de Sitter solution (iTangherlinil Il963l : lGadl2004j ). In 4+1 



for example, the transformation (^j^ now produces a trigonometric expression. 
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3. Extended black hole-de Sitter cosmologies 

For < m^A < i, we construct a few binary black hole spacetimes with different 
separations as shown in Fig. For small separations, the results are very similar to a 

Brill-Lindquist spacetime inside the cosmological event horizon surrounding the binary. For 
large separations, however, the cosmological event horizon splits into two, leading to a binary 
of two universes. Here, the second universe lives in the "tower" of the first, on the sheet 
beyond its cosmological event horizon. 

The dynamical evolution of an extended black hole-de Sitter universe with multiple 
black holes and cosmological event horizons offers a new route to gravitational radiation 
from the early universe. Here, the waves are generated by the dynamics of primordial black 
holes, as well as by multipole moments of cosmological event horizons. The spectrum of 
relic waves may hereby be extended in the infrared, below which what may be expected 
from mergers of primordial black holes alone. It poses a novel computational challenge for 
numerical relativity, i.e., to compute these relic waves at the present epoch from an initial 
distribution of primordial black hole-de Sitter universes at the prior to or at the onset of 
inflation. 



4. Hyperbolicity with curvature-driven lapse 

Spacetime can be describ ed by a foliation in spacelike hyp ersurfaces, given by a 3+1 



decomposition in terms of hij (IThorne. Price fc McDonald! 1 19861 ) 



ds'^ = -N'^dt^ + hij {dx' + I3'dt){dx^ + (3jdt) (12) 

with lapse function and shift functions The and /3* govern the flow of coordinates in 
the foliation of hypersurfaces of constant coordinate time t. In what follows, we shall restrict 
our attention to normal coordinates described by vanishing shift functions, = 0. 

In the Hamiltonian equations of motion with vanishing shift function, the three-metric 
hij evolves according to dthij = —2NKij = —NK^hnj — himNK^. Integration over finite 
time-intervals At gives the product 

KAU) = H,e-^-'"^^*/._(to)e-^"^^*, (13) 

We propose time-evolution with curvature-driven lapse function, i.e., 

DtN = -K, DtKij = -Di,N + }CijN, (14) 
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Fig. 4. — Initial data of binaries constructed by the generalized superposition principle starting 
from zero separation (left top). Shown are isocurves of constant J = \d(p/dx + 1/2 1. For small 
separations, the result is a black hole binary inside a slightly deformed cosmological event horizon 
J = (thick dark line). As the separation increases and becomes large, a binary of two universes 
forms, wherein the second lives on a sheet adjacent to the first and beyond its cosmological event 
horizon corresonding to a dashed sheet in Fig. ([3]) {right below). There is an infinity of black hole 
and cosmological event horizons at exponentially increasing coordinate distances and within the 
two black holes (not shown). 
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where Dij = DiDj — Rij and JC = KKij — 2Kl^Kjm- The gauge condition DtN = —K 
is curvature- driven and is distinct from the product of curvatu re and lapse function in the 
har nionic shcing c ondition dfN = —N^K in Eqs. (69)-(77) of lAbrahams et al.l (Il997l ): see 
also iBrownl ( 120081 ) for a recent review. 

Theorem. The 3+1 Hamiltonian evolution equations are hyperbolic with curvature- 
driven lapse function dfN = —K. 

Proof. It suffices to consider the problem of small amplitude wave-motion about fiat 
spacetime, e.g., the asymptotically flat region with hij = 6ij with = ±1 at large distances. 
Here, we have 



dtN = -K, dlhij = -2Rij + 2DiDjN, d^K = AK. 



We recall that flWaldl 119841 ) 

R 



(15) 



(16) 



where hij = 6hij — j6ij6h, where 6h = h^^dhij refers to the trace of the metric perturbations. 

Small amplitude harmonic perturbations about the fiat metric are given by 5hij ~ 
hije~'^'^^e^^^^\ Conservation of momentum, D^Kij = DjK, the lapse condition —iuN = —K, 
6hij = —2ioj~^Kij and k''Kij = kjK, give rise to 

did^hej — )■ kik'^hej — ^kikj6h = iu}^^{—2kik^Kej + kikjK) = —iuj~^kikjK. (17) 
We then have 

Rij — didjN = -k'^hij — iu'^kikjK + iuj~^kikjK = -k'^hij, (18) 
whereby d'^hij = -2Rij + 2DiDjN gives rise to the dispersion relation 



(19) 



It follows that all small amplitude metric perturbations propagate along the light cone, which 
completes the proof. □ 

Clearly, the system (fTSjl is asymptotically stable, as metric perturbations become small 
at arbitrarily large distances. Asymptotic wave-motion is commonly studied in the so- 
called transverse traceless gauge, or harmonic coordinates-neither of these two coordinate 
conditions are used here. 
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Numerical integration of ( !T3llT^ on the donut (A, 6, (p) can be pursued using a conformal 
decomposition to bring out invariance of the aforementioned Cij with respect to scahng in 
(j). To this end, we may set (p = e^, whereby 



(20) 



It follows that dtKij = Fij with Fij = -N{2K^Kj^ - KK^) + N{Rij + dj) - DiDjN + 
^ijdkN, where fij = 2(j)-^g''^{giedj(f) + gjedicj) - gijdecj)) = 2g^''{giedj'q + gjedif] - gijdev)- 



A closed system is obtained by choosing an equation for the conformal factor. A common 
choice is = h^^^"^ in terms of the determinant h of the three metric hij. For outgoing 
radiation, note that = 1 up to including first order in the wave-amplitude, whereby = 1 
tracks the merger phase of black hole coalescence (here with overlapping horizon surfaces) 
and ringdown in collapse to a single black hole. A slight variation is to insist g^^dtgij = 0. 
The complete hyperbolic system for numerical integration using normal coordinates (/3j = 0) 
hereby becomes 



dtN 
dtV 

dt9ij 



-K, 
-INK, 



-2N ^K,, - y,,K,, 



(21) 



where Kij = e '^^Kij. 

The cosmological event horizon A = A2 represents a turning point in 
an extremum of the radius defined by the Neumann condition 



dfj 
dX 



1 

2' 



representing 



(22) 



As such, fl22|) defines the cosmological event horizon {X2{0,ip),9,{p) following superpositions 
f lTOj) . In dynamical evolutions, it can be used to define apparent horizon surfaces. A proper 
condition of lapse function is that it preserves a simple zero across, i.e., the initial condition 



N 



d log(p0^ 
rfA~~ 



(23) 



5. Conclusions 



We have developed an extension of the black hole-de Sitter space and extended the 
superposition principle for the Hamiltonian energy constraint to (anti-)de Sitter cosmologies 
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with cosmological constant A in terms of an eigenvalue problem for the Laplace operator on 
a metric with constant curvature. 

A positive cosmological constant A is of increasing relevance to our study of cosmological 
spacetimes in view of its role in the early universe in view of observational evidence for a 
small positive value at the present epoch. 

For one black hole, the conformal scale (f) is periodic, representing the extension of the 
black hole-de Sitter space to an infinite tower of universes. The result may become periodic, 
upon identifying any pair of the these event horizons, e.g., Hq with Hi in Fig. 3. Superposi- 
tion of two eigenf unctions produces a binary black hole-de Sitter universe with one (adjacent) 
cosmological event horizon when the separation is small relative to the cosmological scale 
1 / \fK. However, it produces a binary of universes when their separation is large described 
by two adjacent cosmological event horizons. 

The results demonstrate extended topologies of black hole-de Sitter cosmologies, wherein 
the metric is everywhere smooth, as the physical singularities inherent to black holes are 
moved away into the complex plane and space is given a finite physical extent. By virtue 
of a finite evaporation time of black holes by Hawking radiation, there is a strong cosmic 
censorship conjecture by which singularities are protected from direct observation, regardless 
whether they are located in real or complex coordinate space. 

The dynamical interaction of a binary of two universes (j4]) is conceivably of interest 
to the early universe and the ensuing generation of primordial gravitational waves from an 
initial distribution of primordial black hole-de Sitter universes, when A was large at the 
onset of inflation. An extension of our approach to extra dimensions is readily given, which 
is conceivably of interest to dynamical evolution of our four-dimensional spacetime within 
a spacetime of extra dimensions. For computational purposes, we give a new hyperbolic 
formulation of the equations of motion based on a curvature-driven lapse function. 
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